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We describe a method for creating a three-dimensional analogue to Rashba spin-orbit coupling in 
systems of ultracold atoms. This laser induced coupling uses Raman transitions to link four internal 
atomic states with a tetrahedral geometry, and gives rise to a Dirac point that is robust against 
environmental perturbations. We present an exact result showing that such a spin-orbit coupling in 
a fermionic system always gives rise to a molecular bound state. 



Recent experiments with synthetic gauge fields open 
the door to explore spin-orbit coupling and non-Abelian 
gauge fields in atomic systems [TmU] . Non-Abelian gauge 
fields provide rich ground state physics in bosonic sys- 
tems |11H13) . and enhance bound state formation in at- 
tractive fermion systems |14H17) . However, apart from 
engineering cold-atom analogues to known Hamiltonians, 
with suitable choices of laser fields, cold atoms can be 
made to behave in ways that have no known analogue in 
solid state systems 18 . 

In this Letter, we propose a method for synthesizing a 
3D extension to Rashba spin-orbit coupling for ultracold 
atoms which we call Weyl spin-orbit coupling in anal- 
ogy to Weyl fermions [19 . The resulting rotationally 
symmetric dispersion has an infinite ground state degen- 
eracy which covers a sphere, analogous to the Rashba 
Hamiltonian's circular ground state. 

We examine the consequences of such three- 
dimensional spin-orbit couplings: for example, the spec- 
trum for the spherical case has a protected Dirac point 
that cannot be removed by any homogeneous Zeeman 
field. Finally, we consider the addition of interactions and 
exactly show that the 3D spin-orbit coupling (3DS0C) 
strongly enhances the binding energy of two fermions. 

We produce the 3DS0C with a 4-level atom with states 
|1), |2), |3) and |4) optically coupled with a Hamiltonian 

H,i^J2^,k\j){k\ (1) 

with fljk — Q^^'^^ exp[i(kjfc-r-|-0jfc)] as shown in Fig.jlja). 
Here k^j, = — is the momentum transferred by 
the laser, (j)jk is the phase of the coupling, the coupling 
strength is chosen to be rt^^'> for k = j -\- 1, fl^'^^ 

for k — j + 2, and the indexes are taken modulo 4. This 
coupling connects the states in a loop topology "201 , with 
additional next-nearest-neighbor couplings, denoted by 
ri^^). In the maximally symmetric case 
this coupling is geometrically equivalent to a tetrahedron. 
We choose the momentum vectors 

Kj = K±{ex cos /3j — By sin /3j) — K||(— l)"'e2. (2) 



(a) Coupling Scheme 
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FIG. 1. Optical configuration for the 3DS0C. (a) Four states 
are coupled using two-photon optical transitions, (b) The 
momentum-space displacement vectors point to the four 
vertices of a tetrahedron, (c) The fluxes through the surface 
are chosen to be "l>i —71/2 modulo 2n. 



For the remainder of the Letter we will assume the max- 
imally symmetric case with fcy = k±. The vectors Kj 
point from the center to the vertices of a tetrahedron 
[Fig.Hb)]. 

The spatial dependence in the phase term can be elim- 
inated by a state dependent boost \j) — )■ e'^^ '"|j). In the 



2 



boosted basis, the full Hamiltonian is 
(P-K,)2 



(a) States 



2m 



with atom-laser coupling 



(3) 



(4) 
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+ f^(2)^e.fe,+.(|j- + 2)(j|+H.C.) (5) 

The six phases (f>jk are not independent. Only the 
fluxes through the effective surfaces of the tetrahedral 
coupling 



(6) 



are relevant. We choose these fluxes such that $i ~ Tr/2 
modulo 2tt; the sum 

Ei^l = 0, SO of the SIX 0.,t, 

only three are necessary to parametrize the system. We 
use the additional freedom in the choice of phase to elu- 
cidate a symmetry of the problem; in what follows, we 



chose 



7r/4, and 



(j — 1)77 without loss of 



generality. The form of the effective Hamiltonian will be 
the same with any choice of phases for which the fluxes 
through each surface satisfy $i — 7r/2mod27r. 

We diagonalize the atom-laser Hamiltonian through 
a two step process. With the above choice of phases, 
the atom-laser Hamiltonian has a symmetry under the 
transformation |j) — >■ \j + 2). This additional symme- 
try allows us to dimerize states coupled by the next- 
nearest-neighbor couplings, fi^^^ with the transformation 
|a±) = ± |3)) and |&±) = ^(|2) ± |4)). In the 

dimer basis, the atom-laser Hamiltonian is 



Ha 





V 



/2nw 



i\/2r2(i) 






-^V2^^(l) 



(7) 



with respect to the basis {|a+), |5+), |a_)}. The 
subsequent unitary transformation U independently di- 
agonalizes the two blocks, which have the same spec- 
trum, E ^ zt^ 



'2 [nW] + [17(2)] . The twofold degener- 
ate ground states comprise our pseudospin for the spin- 
orbit coupling. 

The spin-orbit coupling arises from projecting J2j Q ' 
Kj\j){j\ into the low energy subspace to give q- A, where 
the 3D vector potential 



{(7xex + o-yBy) + sin^Kiicr^e^ 



(8) 



nontrivially includes all three components of the Pauli 
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(c) Laser Geometry 
Li4 



BEE T ^ 

^24^34 — ^ 



-t* ■<»- 



12-t^l3 



e_ 



FIG. 2. Four hyperfine states |_F, wif) of *^Rb are coupled 
using nine lasers. The quantization axis is set by a Zeeman 
field along the S-axis. The coupHngs are produced in pairs, 
(a) The four states in the tetrahedral coupling are mapped 
to physical states according to |1) — j2,0), |2) — 
jS) = |1, 0), |1) = |2, +1). (b) The frequencies of the three sets 
of lasers are given by {uJa,ijJa + Si2,uJa + S34,} (dashed blue), 
{ujb,ujb + Si3,ijJb + $24} (dotted black), and {cJc, i^c + 5i4,ujc + 
S23} (solid red), where Sij — toi—ojj is the frequency difference 
between the states \i) and \ j) in the rotating frame, (c) The 
geometry of the nine laser beams. — ki2 = — ki3 — k23 = 
k34 = fc^e-, ka = k23 = — ki4 = kLZ and k;, = — kc = fciC-. 
The unit vectors e± — ±^ (izby). For a complete description 
of the laser parameters, see the Supplemental Material. [21] 



17(2) /217(^) and , the vector potential can give both 

symmetric and asymmetric spin-orbit couplings. The 
3DS0C is fully isotropic when kj_/k\\ = n^^^ /2n^^\ with 
a Hamiltonian 



Hn 



P 

2m 



vcr ■ p. 



(9) 



matrices {(Jx,o-y,az) 



By changing tan 9 



The spin-orbit coupling is characterized by the velocity 
V = Keff/rn, where Ate// = k± cos(20)/2. 

The 3DS0C can be implemented in ^^Rb using two- 
photon transitions. A possible implementation is given 
in Fig. [2] Nine laser beams with wavelength A are used 
to couple states within the F — 1 and F = 2 hyperfine 
manifolds. A Zeeman field of i? = 200mT sets the quanti- 
zation axis along the z direction. The remaining hyper- 
fine transitions are isolated with a 6.8 GHz microwave 
field. Each pair of nonadjacent couplings is induced with 
three laser beams. For example, the optical couplings 
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(a) Ground state manifold 

E=-mv"/2 




(b) Energy shells 




(c) Cross-section 

E(p) 




FIG. 3. Energy spectrum of the isotropic 3DS0C. (a) Energy 
spectrum of spherical ground state manifold |p| = mv with 
E = — mii^/2. Near the ground state manifold the disper- 
sion E{Sp) is approximately parabolic in a small deviation, 
Sp, from the momentum of the ground state, (b) Shells of 
constant energy, E„ — —(1 — 0.1n)mu^/2 with n = 0, 1,2,3. 
These low energy shells have smaller surface area until p = 0. 
(c) A one dimensional cut of the energy spectrum. The full 
spectrum is generated by rotating this spectrum along three 
ajces. The blue(red) band corresponds to states with momen- 
tum aligned (antialigned) with spin. 



rii2 and r234 are produced using the beams La, L12 and 
L34. These beams have the respective frequencies uja, 
cJq + wi — LU2 and -|- wa — W4 , and corresponding polar- 
ization vectors — '^{^ ~ y)j ^12 = ^34 = ^- This pat- 
tern will ensure the system does not undergo unwanted 
optical transitions. The remaining four optical couplings 
are shown in Fig. [2j For complete details see the Supple- 
mental Materials. [21 

The spectrum of Eq. [9] is given by 



Eip) 



2m 



±vJpl 



(10) 



and is shown in Fig |3] At p = 0, the system has a 
protected three-dimensional Dirac point. Since the spin- 
orbit coupling includes all spin matrices, no uniform Zee- 
man field can induce a splitting at the Dirac point. This 
is true even for anisotropic couplings, provided that the 
nature of the 3D spin-orbit term is preserved. 

The spectrum has a ground state manifold on a sphere 
IpI = mv. The dispersion is parabolic only along the ra- 
dial direction, and excitations are energetically free along 



the polar and azimuthal directions. This suggests that 
a 3DS0C will behave as a quasi-one dimensional sys- 
tem with respect to fluctuations around the ground state. 
This will be manifest in the presence of a spherically sym- 
metric trapping potential. Provided the spin-orbit energy 
is sufficiently strong, the low energy spectrum will be de- 
fined by a single radial quantum number, and will be 
degenerate in the angular quantum numbers. [22j 

Similarly, two free fermions with spin-orbit coupling 
will have a binding energy that is enhanced. It is well 
known that in one and two dimensions, an arbitrarily 
weak, attractive, potential has a bound state, whereas 
no bound state is guaranteed to exist in three dimen- 
sions |23j . Furthermore, the binding energy in a ID 
system is algebraic in the potential, while in a 2D sys- 
tem the binding energy is exponentially small. It was 
noticed [Ml , and recently rediscovered [TiHTB] , that 
Rashba fermions have an enhanced tendency for molec- 
ular formation. This can be understood in a manner 
similar to the formation of Cooper pairs, where the en- 
hancement of the density of states near the Fermi surface 
reduces the effective dimension of the system from three 
to two. 

The situation is analogous for 3DS0C fermions, where 
the energetically free excitations along the polar and az- 
imuthal direction enhances the density of states, and 
effectively reduces the bound state problem to one di- 
mension. We now present the summary of an exact 
calculation demonstrating the binding energy of two 3D 
spin-orbit coupled Fermions is enhanced. For the com- 
plete calculation see the Supplemental Materials [35] . To 
search for bound states of two Fermions with 3D spin- 
orbit coupling, we solve the two-particle Schrodinger's 
equation 

[i/,o(ki) ® i + i $5 i7,o(k2)]|*) + ^12!*) = A|*), (11) 

where V12 is the two-particle interaction potential, and 
the tensor product implies an operator on the left op- 
erates on particle 1, and the operator on the right op- 
erates on particle 2. We assume the interaction to be 
purely local s-wave. This equation can be expressed in 
self-consistent form as 



|*> =GT/i2|*) 
where the Green's function is defined as 



(12) 



G = [{Hso - A/2) ®i + i®{H,o- A/2)]~i. (13) 

The ground state of a spin-orbit coupled atom will have 
energy Eso — —mv'^/2. We therefore search for solutions 
of (171 with energy E < 2Eso = —mv"^. We define the 
binding energy as e = —rav^ — A > 0. 

The s-wave character of V12 will project the ground 
state into the singlet channel. Upon integration over the 



4 



relative momentum coordinate k = ki 
press the self-consistency equation as 



we can ex- 



d^k 

(27r)3 



{^s\G{-k,QMs)^l, 



(14) 



where \ips) is the projection of the exact ground state 
into the singlet channel, Q = ^(ki + k2) is the center 
of mass momentum, and vq is the characteristic scale of 
the interaction. The bound state |^'b(Q)) will appear as 

a solution to this equation with energy A < — 2 ^^^^ , 

which is twice the ground state of of a single spin-orbit 
coupled fermion. We find that the binding energy of two 
fermions is given by 



£ = A- 



mv = —mv 



(15) 



at Q = to lowest order in e/mv^. Thus, there exists 
a negative energy two-fermion bound state with energy 
that is algebraic in the interaction strength, consistent 
with the mean- field results obtained previously [14] . This 
is in contrast to a system of two fermions without spin- 
orbit coupling where there is no bound state in three 
dimensions for weak attraction. The existence of this 
bound state results from the enhanced density of states 
near the 2D surface defined by |p| = mv. The addi- 
tional states increase the effect of quantum fluctuations 
and provide an effective dimensional reduction of the sys- 
tem by 2 to D = 3 - 2 = 1. 

For nonzero center of mass momentum Q, the binding 
energy of the system becomes Q dependent. For Q <C 
mv the binding energy is e(Q) = —mv'^{mvao)^ ~ 3^ + 
0{Q)'^, which is the dispersion for the center of mass of 
the free particle. At large momenta near Q ^ mv, the 
self-consistency equation no longer has a solution, and no 
bound state is possible. We note that the binding energy 
is a monotonically decreasing function of A, and thus 
there is exactly one bound state for sufficiently small Q. 
This single bound state results from the delta function 
interaction. Additional bound states may be possible 
for other interactions, such as a square well or a p-wave 
interaction. 

In conclusion, we proposed a scheme to produce a 
3DS0C using two-photon transitions to couple four 
atomic levels in a tetrahedral topology. In the limit 
of large optical power, the ground state is defined by 
a sphere. The origin of this coupling can be viewed as an 
approximation of the desired spherical momentum-space 
ground state manifold by a tetrahedron, whose spherical 
symmetry is restored in the infinite coupling limit. Such 
a coupling could give rise to interesting many-body sys- 
tems, such as a Bose liquid [27, or with the addition of 
suitable band gaps, Weyl fermions. Finally, we present an 
exact solution for the bound state energy of two fermions 
with a spherical 3D spin-orbit coupling. Such a bound 
state is found to always exist for sufficiently small center 



of mass momentum, and the energy of the bound state 
is algebraic in the interaction strength. 
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Note: After the completion of this work, we were made 
aware of a similar calculation on the bound states of 
3DS0C coupled fermions but not the possible 

origin of such coupling, as well as a possible implemen- 
tation of 3DS0C on a lattice. [30] 
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SUPPLEMENTAL MATERIALS: 
BOUND STATES 



FERMIONIC 



To search for bound states of two Fermions with 
3D spin-orbit couphng, we solve the two-particle 
Schrodinger's equation 

[i/,o(ki)®i-|-i®i/,„(k2)]|*) -ffial*) = A|*), (16) 

where V12 is the two particle interaction potential, and 
the tensor product implies an operator the left operates 
on particle 1, and the operator on the right operates on 
particle 2. We assume the interaction to be purely local 
s-wave. This equation can be expressed in self-consistent 
form as 

\^)^GVi2\^) (17) 

where the Green's function is defined as 

G = [{H,o " A/2) ® i + i ® (i7,„ ~ A/2)]-\ (18) 

The ground state of a spin-orbit coupled atom will have 
energy Ego — — mu^/2. We therefore search for solutions 
of (171 with energy E < 2Eso — —mv^ . We define the 
binding energy as e = —mv^ — A > 0. 

The Green's function can be calculated by applying 
the unitary matrix U — Ui®U2, where 



Uj = exp 



(19) 



is the unitary matrix that rotates from the original 
spin basis to the pseudo-spin basis defined by |a/3), 
with a,/3 = ±1, where a particle of spin |±) has en- 
ergy in the E = Tp^ /2m ± t)|p| band. The vector 
(— sin0j, cos(?!)j,0), is perpendicular to both kj = 
cos (j) j, cos 9 jSmOj^smOj) and e^. The unitary 



3 

kj (cos 

matrix Uj transforms the Green's function to 



(20) 



al3 



where dap — [s + v{aki + /?fc2)] 

We now assume the interaction potential is a 
short range, s-wave interaction, V12 = —voS{ri — 
T^2)'Ps, where Vs is a projector into the singlet state 
\i^s) = |(in>-|;t>)- Using (|l7|, we first ap- 
ply the interaction potential to me state j^I') to 

get V^2{\A^Bm = / (0(Vs,k|vI/B(Q))|V.), where 



(k|^'B(Q)) is the momentum-space wavefunction of the 
relative coordinate k — (ki — k2)/2. The center-of-mass 
momentum Q ~ ki + k2 commutes with the Hamil- 
tonain, and is thus a good quantum number which la- 
bels the bound state |\E'b(Q)). We can therefore ex- 
press |*b(Q)) = A^(Q)G(k,Q)|V^), where iV(Q) 
/(0(V'.,k|^PB(Q)). 



We then find the full wavefunction by applying the 
Green's function to the state |5'(Q)) to get 



IVb(Q)) = iV(Q) "^^^ 



(21) 



where q;,/3 = ±1, and the spinors Xap are most easily 
expressed in the coordinates ki , k2 as 



Xafi 



/ il3 + COS 9i) sin 626-"^'^ + {a + cos 6*2) sin die-"^^ 
(a -I- cos 0i){-f3 + cos 62) - sin 0i sin 02e''-'t'^-'f'^'> 
-(-a + cos6'i)(/3 -h COS6I2) + sin6li sin6'2e-^('^^-'^i 
V (-/3 4-0086(2) sin 6*1 e*"^! (-a + cos6li) sin6'2e**2 

(22) 

in the basis of (|tt), lU), l^t), |U))^- The form 
an orthogonal basis for each ki , k2 , but are not nor- 
malized, since xLp'Xap = 
8 {1 + a/3 [cos 6*1 cos 62 + cos( 

In this notation, the normalization is -/V(Q) — 
The full wavefunction allows us to calculate the bind- 



6aa'Sl3l3'Cal3, and Cq^ 

i>i - 02) sin 01 sin 02]}- 



ing energy as follows. We substitute (21 1 into (17 1, and 



then integrate over the relative momentum coordinate. 



Due to symmetry, the triplet components of (21) will 



vanish under integration over the relative momentum, 
/ dik ^klV'B(Q)) = / 707'.(k|^s(Q)). Therefore, left 



multiplying (17) by {i<-\Vs, and integrating over k allows 



us to express the self-consistency equation as 
d=^k 



i2ny 



r(V's|G(k,Q)|V's) = 1. 



(23) 



The bound state |^_b(Q)) will appear as a solution to 

this equation with energy A < — 2 ^^^^^ , which is twice 

the ground state of of a single spin-orbit coupled fermion. 
The integrals over the four modes dap have a linear ultra- 
violet divergence. To regularize the integrals we replace 
the interaction with a renormalized scattering length 



1 



1 



1 



(24) 



where l/woo = J2a0 1 



^0 J (27r)3 (fc-|-Qm«)2 + (fe+^mi,)2 ■ ^ 

form of this regularizer is necessary to cure the lin- 
ear divergence. A regularization scheme of 1/woo 
I m/k'^ will reduce the divergence from linear to 
logarithmic. Calculating the integrals using this regu- 
larization scheme above can be performed exactly, giving 
a binding energy 



A 



-mv (mvao) 



(25) 



at Q = to lowest order in e/mv^ . 
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SUPPLEMENTAL MATERIALS: 
IMPLEMENTATION 

The 3DS0C can be implemented in *^Rb using two 
photon transitions. A possible implementation is given 
in Fig. [2] Nine laser beams with wavelength A are used to 
couple states within the F — 1 and F — 2 hyperfine man- 
ifolds. A Zeeman field B — 200mT sets the quantiza- 
tion axis along the z direction. The remaining hyperfine 
transitions are isolated with a 6.8 GHz microwave field. 
Each pair of non-adjacent couplings is induced with three 
laser beams. All lasers have frequencies tuned between 
the Di and D2 transitions. To ensure that all unwanted 
couplings are off resonance, we chose the base frequen- 
cies much larger than the splitting between the hyperfine 
levels Wq, ojb, ojc ^ Sij . 

The six optical couplings are induced in non-adjacent 
pairs on the state-linkage diagram. We list the properties 
of the three pairs of couplings independently. In what 
follows we define the vectors e± = {x ± y) and hk^ 
is the recoil momentum of the laser. 



Couplings r2i2 and ^34 

The two couplings fli2 and fl^^ will be induced with 
three lasers, denoted La, L12 and L34. The beam La will 
be shared in the two-photon couplings. The frequencies 
of the three beams will be chosen such that 



(a) States 



1^34 



(5l2 



(26) 
(27) 



where Ua is the frequency of the laser La and 5ij , is the 
frequency splitting between the states \i) and \i) . With 
such a configuration, the transition |1) o |2) and |3) O 
1 4) will be on resonance, while all other dipole allowed 
transitions will be off resonance. The wavevectors of the 
lasers are given by 

K.a = kLZ (28) 
«12 = -kLe- (29) 
«34 = fcLe_. (30) 

We can check that the effective momentum transfer of 
the couplings are given by 



'■12 



Ki K2 



\/2z) 



and 



^■34 



fcL(-e- +i) 

K3 - K4. 



V2i 



(31) 
(32) 

(33) 

(34) 
(35) 
(36) 



(b) Laser Couplings 
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(c) Laser Geometry 
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FIG. 4. Four hyperfine states \F,mF) of *''Rb are coupled 
using nine lasers. The quantization axis is set by a Zeeman 
field along the S-axis. The coupHngs are produced in pairs, 
(a) The four states in the tetrahedral coupling are mapped 
to physical states according to |1) = |2,0), |2) = 
|3) = |1, 0), |1) = |2, +1). (b) The frequencies of the three sets 
of lasers are given by {tJa,'^a -I- <5i2,i^a -I- 52.4,} (dashed blue), 
{tJt, ojb + (5i3, cjb + 524} (dotted black), and {cJc, i^c + (5i4, ti^c + 
52-i} (solid red), where 5ij — toi—ujj is the frequency difference 
between the states \i) and \ j) in the rotating frame, (c) The 
geometry of the nine laser beams. — ki2 = — kia — k23 = 
k34 = kLe-, ka = k23 = — ki4 = kLZ and k(, = — k^ = kLC-. 
The unit vectors e± — ±^(i: ± ?/). 



Finally, the polarizations of the lasers will be chosen such 
that the beams L12 and L34 are linearly polarized along 
the z direction, while La is (t+ polarized. 



Couplings Qi4 and ^23 

The couplings f2i4 and ^23 will similarly be induced 
with three lasers, denoted Lf,, L14 and L23- The beam 
Lh will be shared in the two-photon couplings. The fre- 
quencies of the three beams will be chosen such that 



Wl4 
W23 



Si4 



(37) 
(38) 



where Wf, is the frequency of the laser Lf,. This choice 
of frequencies will isolate the transitions |1) O |4) and 
|2) o |3) in a manner similar to above couplings. The 
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wavevectors of the lasers are given by 



Kb 
1^23 



kLe+ 
khZ. 



(39) 
(40) 
(41) 



We can check that the effective momentum transfer of 
the couplings are given by 



. kL{e+ + z 
kL^ 

Ki - K4 



y- 



and 



k23 = kL{e+ - z) 



kL . 
= Ko - K, 



(42) 
(43) 
(44) 

(45) 
(46) 
(47) 



Finally, the polarizations of the lasers will be chosen such 
that the beams Li, is linearly polarized along the z direc- 
tion, while 1/14 and L23 are linearly polarized along the 
X axis. 



Couplings f2i3 and Q,2i 



The last pair of couplings f2i3 and f224 will similarly 
be induced with three lasers, denoted Lc, iis and i24- 
The beam will be shared in the two-photon couplings. 
The frequencies of the three beams will be chosen such 
that 



Wl3 = Wc + <5l3 
W24 = Wc + ^24, 



(48) 
(49) 



where lJc is the frequency of the laser Lc- This choice 
of frequencies will isolate the transitions |1) O |3) and 
1 2) -H- |4) in a manner similar to above couplings. The 
wavevectors of the lasers are given by 



«24 



-kLC- 



kre- 



(50) 
(51) 
(52) 



We can check that the effective momentum transfer of 
the couplings are given by 



tl3 



kL{-e+ + e-) 

K1-K3 



(53) 

(54) 
(55) 



and 



ko4 = 



(56) 
(57) 



kL{-e+ - e_) 

= K2 - K4. (58) 
Finally, the polarizations of the lasers will be chosen such 
that the beams L12 and L34 are linearly polarized along 
the e_|_ direction, while is linearly polarized along e_ . 



Amplitude and Phase 

In each pair of transitions, the three lasers provide a 

sufficient number of both amplitude and phase degrees 
of freedom to chose the the values of the couplings as 
desired in the main text. 



